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We present two issues here: (i) that in situation when total energy is kept constant re- 
cently proposed fluctuations of volume ensemble is equivalent to the approach using Tsallis 
statistics with fluctuating temperature and (ii) that the later (in which fluctuations are 
described by the nonextensivity parameter q) leads to the observed experimentally sum 
rule connecting fluctuations of different physical observables. 

1 Introduction 

Statistical modelling represents a standard tool widely used to analyze multiparticle production 
processes pQ. However, this approach does not account for the possible intrinsic nonstatisti- 
cal fluctuations in the hadronizing system which usually result in a characteristic power-like 
behavior of single particle spectra or in the broadening of the corresponding multiplicity dis- 
tributions (and which can signal a possible phase transition(s) [2]). To include such features 
one should base this modelling on the so called Tsallis statistics [3J 21 [SJ (represented by Tsallis 
distribution) which accounts for such situations by introducing, in addition to the temperature 
T, one new parameter, q > 1, directly connected to fluctuations [6l [7] (for q — > 1 one recovers 
the usual Boltzmann-Gibbs distribution): 
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The most recent applications of this approach come from PHENIX Collaboration at RHIC [8] 
and from CMS Collaboration at LHC [5J. One must admit at this point that this approach is 
subjected to a rather hot debate of whether it is consistent with the equilibrium thermodynamics 
or it is only a handy way to phenomenologically description of some intrinsic fluctuations in the 
system under consideration [10] . However, as was recently demonstrated on general grounds in 
[llj , fluctuation phenomena can be incorporated into traditional presentation of thermodynamic 
and Tsallis distribution [3J belongs to the class of general admissible distributions which satisfy 
thermodynamical consistency conditions and which are therefore a natural extension of the usual 
Boltzman-Gibbs canonical distribution. Actually, what was shown in [3] was that starting from 
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some simple diffusion picture of temperature equalization in the nonhomogencous heat bath (in 
which local T fluctuates from point to point around some equilibrium temperature, To) one gets 
evolution of T in the form of Langevin stochastic equation and distribution of 1/T, f(l/T), as 
solution of the corresponding Fokker-Planck equation. It turns out that /(1/T) has form of 
gamma distribution, 

Convoluting exp(— E/T) with such /(1/T) one gets immediately Tsallis distribution, exp 9 (— E/T) 
from Eq. (JTJ [5J. Parameter q, i.e., according to Eq. ((T|) also the temperature fluctuation pat- 
tern, is therefore fully given by the parameters describing this basic diffusion process (cf., [B] 
for details, this was recently generalized to account for the possibility of transferring energy 
from/to heat bath, which appears to be important for AA applications [H[12] and for cosmic 
ray physics [13] ; we shall not discuss this issue here) . This approach has now been successfully 
applied in many circumstances, see [H Q21 El [9] and references therein. 



2 Fluctuations of V or T? 

It must be stressed at this point that the form of /(1/T) as given by Eq. ([2]) is not assumed but 
has been derived from the underlying physical process. We shall now compare this approach 
with that proposed in |14j in which the volume V was assumed to fluctuate in the scale invariant 
way following the observed KNO scaling behavior of the multiplicity distributions, P(N) (15] . 
We shall demonstrate here that when total energy is kept constant, as was assumed in [14] . 
both approaches are equivalent. Let us first notice that for constant total energy, E = const, 
both the volume V and temperature T are related via E ~ VT A , what means that 

Following now [14) the mean multiplicity in the microcanonical ensemble (MCE), N, can be 
written as 

V ( T N 



N = {N) -(Y)W)) = {N)V > (4) 

what means that N fluctuates in the same way y. It is then natural to assume that y follows the 
pattern of fluctuations of N, i.e., KNO limit of the NBD distribution observed in data fitting 
|15j . which is given by gamma function. The power-like form of single particle spectra then 
follow immediately, all apparently without invoking any reference to Tsallis statistics. Notice, 
however, that because of ([3]), 1/T will fluctuate according to the same gamma distribution. So 
we get T fluctuations with the same functional form but now without the physical background 
behind the Eq. ([2} mentioned above. However, we can proceed in reverse order and obtain 
from our T fluctuations introduced in Section [T] fluctuations of V introduced in [14) . In this 
sense both approaches are equivalent with the former being based on some physical processes 
and the second on apparently ad hoc assumption^. 



however, after all, this assumption can have some phenomcnological foundation, not mentioned in [14] . which 
deserves further scrutiny. Namely, one observes experimentally a variation of the emitting radius (evaluated from 
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We close this section with short reminder that temperature fluctuations discussed in Section 
[1] result in automatic broadening of the corresponding multiplicity distributions, P(N), from 
the poissonian form for exponential distributions to the negative binomial (NB) form for Tsallis 
distributions |17j . It is known that whenever we have TV independently produced secondaries 
with energies {J3i=i .,jv} taken from the exponential distribution in Eq. ([I]) and whenever 
Y^Ln Ei < E < X^ilo 1 Ei, then the corresponding multiplicity distribution is poissonian, 
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What was shown in [T7] is that whenever in some process N particles with energies {Ei—i 
are distributed according to the joint TV-particle Tsallis distribution, 
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(for which the corresponding one particle Tsallis distribution function in Eq. ([T]), is marginal 
distribution), then, under the same condition as above, the corresponding multiplicity distri- 
bution is the NB distribution, 
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Notice that in the limiting cases of q — > 1 one has k — > oo and ([7]) becomes a poissonian 
distribution (J5)), whereas for q — > 2 on has k — > 1 and ([7]) becomes a geometrical distribution. 
It is easy to show that for large values of TV and (TV) one obtains from Eq. its scaling form, 
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and 



in which one recognizes a particular expression of Koba-Nielsen-Olesen (KNO) scaling 
which, as discussed before, has been assumed to describe also the volume fluctuations in pH] . 
Here it results from the temperature fluctuations described by the parameter q discussed in 
Section [T] with well defined physical meaning. 



3 Composition of different fluctuations 

Description of fluctuations phenomena by means of the parameter q using Tsallis statistics 
allows for better understanding of interrelations between different fluctuations. From our ex- 
perience with pp collisions |19| we know that one can obtain very good description of the 

the Bose-Einstcin correlation analysis) with the charged multiplicity of the event, see, for example, 1161 . An 
increase of about 10 % of the radius when the multiplicity increases from 10 to 40 charged hadrons in the 
final state was reported. Unfortunately, the quality of data does not allow us to precisely determine the power 
index of the volume dependence. It is also remarkable that both the energy density, = E/V, and particle 
density, pjv = N/V, decrease for large multiplicity events. For N/(N) ~ y one observes Pe/{pe) ~ V~ 4 and 
Pn/{pn) ~ V All these deserves further consideration and should be checked in future LHC experiments, 
especially in ALICE, which is dedicated for heavy ion collision. 

2 It is worth to mention at this point that, as shown in 118] . fluctuations of TV in the poissonian distribution 
taken in the form of ij)(N / < N >), Eq. JSJ, lead to the NB distribution Q. 
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whole range of p T (oc cxp q (~p T /T) with (T T [GcV];q T ) = (0.134; 1.095), (0.135; 1.105) and 
(0.14; 1.11) for energies (in GcV) 200, 540 and 900, respectively. These values should be com- 
pared with the corresponding values of (T = Tl; q = qi) obtained when fitting rapidity distribu- 
tions (oc exp g (-fir coshy/T)) at the same energies: (11.74; 1.2), (20.39; 1.26) and (30.79; 1.29). 
It was noticed there that qi, — 1 has the same energy behavior as 1/k in the NB distribution 
fitting the multiplicity distributions at corresponding energies (q^ — 1 = — 0.104 + 0. 058 In -^/s). 
It means that fluctuations of total energy are in this case driven mainly by fluctuations in 
the longitudinal phase space. Explanations proposed in |19j was following. Noticing that 
q — 1 = a 2 (T)/T 2 (i.e., it is given by fluctuations of total temperature T) and assuming that 
<7 2 (T) — (t 2 (Tl) + <t 2 (Tt ) , one can estimate that resulting values of q should not be too different 
from 

q L T 2 + q T T 2 T 2 +T 2 TlMt 

1 = f2 f2 + 1 ~ 9L- (9) 

It turns out that situation is completely reversed in the case of nuclear collisions, which we 
shall discuss now{f|, cf. Fig. [1] Left panel shows q obtained from different sources as function of 
centrality represented by number of participants, Np. The one obtained from P(N) follows 

behavior (a = 0.98) [T^]. Whereas for small centralities it approaches situation encountered in 
pp collisions (where it was practically equal to q = qh obtained from rapidity distributions as 
mentioned above), the more central is event the smaller is q — 1, i.e., the nearer to poissonian 
is the corresponding P(N). Notice that both, and qr (obtained from px distributions arc 
now greater than q and have (approximately) visible similar dependence on Np, however now 
qh < <?Ti again opposite to what was seen in pp [19) (for comparison qT obtained by |21j using 
RHIC data Au + Au collisions at 200 GeV [33] are shown here as well) . Right panel shows the 
same quantities but now as function of energy for the most central Pb + Pb collisions [20] . In 
both cases we take from [20] distributions of rapidity, dN/dy, and in \ij<, dN/d^r, and from 
them deduces the corresponding q. 

The natural question is, what causes such different behavior of parameter q in this case. 
The answer we propose is the following. When extracting values of parameter q from the 
rapidity distributions tacit assumption was made that {It in E = fir cosh y remains constant 
(i.e., it does not fluctuate). What would happen if this assumption was false? Notice that in 
the exp ? (-E/T) = exp q [— (px/T) coshy] = cxp 9 (— z coshy). It means that fits to rapidity 
distributions provide us, in fact fluctuations not so much of partition temperature T but rather 
of the variable z = [i^jT. This in turn can be written approximately as: 

tt i \ 1 r , ,^^T? Var(T) 

Var(z) ~ -jrjyVar ( M t) + -^p j^^- (11) 

Because (z) ~ (/x T )/(T) and Var(l/T)/(l/T) 2 ~ Var(T)/(T) 2 and because Var{z)/(z) 2 = 
Var {(i T ) /(^t) 2 + Var(T)/(T) 2 one can write that 

_ , dg Var(T) = Var(z) _ Var (// T ) 
q (T) 2 (z) 2 ( MT )2 ■ ^ ' 



3 We use for this purpose NA49 data on Pb + Pb collisions 1201 because, at the moment, only this experiment 
measures at the same time (at least for the most central collisions) multiplicity distributions, P(N), and dis- 
tributions in rapidity y, transverse momenta, pt, and transverse masses, = y m? + •pip , which is crucial for 
our further considerations here. 
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Figure 1: Left panel: q from P(N) are our results obtained in [T2] from Var(N)/(N), q from 
RHIC are taken from compilation [21] and are based on data for dN/dpT from [22] analyzed 
there, whereas q from NA49 are obtained by us for this presentation using data on dN/dpx from 
the first work in [2D]. Right panel: all results were obtained for the sake of this presentation 
using distributions provided by [20] . i.e., respectively, dN/dpT, dN/dy and dN/dpT- The errors 
are similar to those presented as example for q obtained from dN/dy. Open symbols correspond 
to uncorrected values of q, full symbols to values corrected by means of the procedure proposed 
in the text. 



This sum rule is our main result and its action is presented in the right panel of Fig. [T] 
It connects total q, which can be obtained from the analysis of the NB form of the measured 
multiplicity distributions, P(N), with qi, — 1 = Var(z)/{z) 2 , obtained from fitting rapidity 
distributions and Var (pr) / (a*t) 2 obtained from data on transverse mass distributions. When 
extracting q from distributions of dN/ d^ir we proceed in analogously way with z being in this 
case equal to z = coshy/T. 



4 Summary 

To summarize: we have demonstrated that for constant total energy fluctuations of T introduced 
by us some time ago [6] [4] are equivalent to fluctuations of V proposed recently [14] and that, 
at the moment, the former have advantage of being backed by a plausible physical arguments 
[DJ. Moreover, due to relation ([3]) valid for constant total energy, the inverse temperature 1/T 
and I/ 1 / 4 fluctuate in the same way, according to gamma distribution and such fluctuations 
lead to the Tsallis form of the respective distributions for energy spectra. 

However, the main results presented here is the sum rule formula, Eq. (|12j) . connecting q 
obtained from analysis of different distributions which are obtained in the same experiment. 
This allows us to understand why in AA collisions fluctuations observed in multiplicity distri- 
butions are much smaller than the corresponding ones seen in the rapidity distribution or in 
distribution of transverse momenta (i.e., why the corresponding q parameters evaluated from 
distributions of different observables are different). This issue should be checked further when 
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completely sets of data would become available from the experiments at LHC (especially from 
ALICE). 
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